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An explanation of the difference of the charge radius of the proton as determined from the Lamb 
shift in electronic hydrogen and from elastic electron scattering off the proton on the one side and 
the recent high precision determination with muonic hydrogen on the other side is presented, ft is 
shown that the modification of the 2Si/2 and 2P 3 / 2 wave functions by the "Uehling potential" yields 
a correction to the theoretical Lamb shift of 5(A_ELamb) = 0.302 meV which has to be compared to 
5(Ai?Lamb) = 0.322(46) meV equivalent to the stated radius difference. The explanation is based on 
the realization that the bound state wave functions modified by the external "Uehling potential" 
have to be propagated by the vacuum polarization propagator in order to give the correct leading 
order Lamb shift. It is argued that a conflicting relativistic calculation neglects this propagation 
aspect. The explanation demonstrates that the Lamb shift is dynamically induced through the QED 
vacuum polarization and is not only the result of a static external "Uehling potential" probed by a 
test charge. 

PACS numbers: 14.20.Dh, 31.30.jr, 36.10.Ee , 25.30.Bf 



INTRODUCTION 

The seven standard deviations difference of the root- 
mean-square (rms) charge radius of the proton r p = 
( r p) 1 ^ 2 derived from electron scattering and from the 
Lamb shift of muonic hydrogen has caused a consider- 
able worry in the physics community. Since the Lamb 
shift is a corner stone of the tests of QED this difference 
requires indeed a convincing clarification. In this letter a 
proposal for an explanation in the framework of standard 
QED is given. 

The rms radius derived from the Lamb shift in muonic 
hydrogen is r p — 0.84184(67) fm [1] which has to be com- 
pared to the CODATA value of 0.8768(69) fm [2] mean- 
ing a five standard deviations difference. The CODATA 
value is derived from a measurement of the Lamb shift 
in electronic hydrogen [3] and from electron scattering 
experiments |4J. Almost at the same time as the muonic 
experiment a new independent determination with elas- 
tic electron scattering has been published and yielded 
0.879(8) fm |5j resulting in a weighted average of all "elec- 
tronic experiments" of 0.878(5), i.e. the mentioned seven 
standard deviations difference. 

In order to transform this radius difference into an en- 
ergy deviation of the muonic Lamb shift we repeat the 
key formula for the Lamb shift pQ: 

AE hZb = (209.9779(49) - 5.2262 ^/fm 2 + 

0.0347 r3/fm 3 )meV (1) 

The experiment yielded A^Xb = 206.2949(32) meV 
from which the radius was determined. The term with r p 
is an approximation to the 3rd Zemach moment (r 3 )^) 
depending on the charge form factor of the proton. In- 
serting the differing radii and their error bars given above 



one gets for the deviation of the Lamb shift: 

6(AE Lamh ) = 0.322(46) meV (2) 

where the error is dominated by the electronic experi- 
ments. 

However, a recent non-perturbative relativistic cal- 
culation of the theoretical Lamb shift [6] yields some- 
what different constants for the formula in eq. ([l}. From 
this follows r p = 0.83340(67) fm resulting in a deviation 
(5(Ai?Lamb) = 0.403(51) meV. Considering the perfection 
of this calculation one is puzzled about possibilities of an 
explanation. 

In view of the excellent accuracy of the muonic ex- 
periment and the very good agreement of the about half 
dozen electronic experiments, it is highly unlikely that 
the reason for the deviation is due to a problem on the 
experimental side. Therefore, one has to see what possi- 
bilities are left in the analysis of the data and in the frame 
work of the QED description. At first one has to realize 
that the QED calculations of the muonic Lamb shift have 
been scrutinized again and again over the years. Of the 
many publications we mention some pertinent summaries 
[THTTj. On the other hand, there is also no room for a 
modification of the radiative corrections to the electron 
scattering cross sections. 

In this paper it is shown that the wave functions deter- 
mined by Carroll et al. [6] have been used incompletely 
as input in the QED description of the Lamb shift. If 
correctly used the radius puzzle disappears. 

PREVIOUS PROPOSALS FOR EXPLAINING 
THE DIFFERENCE 

Of the many proposals for explanations we list some 
prominent and typical. 
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• 3rd Zemach moment 

One can construct an electric form factor with an 
evanescent charge cloud extending to very large 
radii but maintaining the radius of the electric ex- 
periments. This increases the 3rd Zemach moment 
( r p)(2) so much that the muonic radius determined 
from cq.Q agrees with the electric one [T5Ml4j . 
However, it has been shown that such a conjecture 
is in disagreement with the measured form factors 

QMS]. 

• Off-shell form factors 

Miller et al. (TU] have proposed off-shell influences 
of the form factors in the elastic box diagram con- 
tributing to the muonic hydrogen Lamb shift. This 
idea was refuted by Carlson and Vanderhaeghen 
who showed that such contributions are two orders 
of magnitude smaller than thought [2"0] , 

• Form factor extrapolation problems 

The radius is derived from the electric and mag- 
netic form factors Gem according to (r 2 ) 1 / 2 = 
(6/G EM {Q 2 )){dG EM (Q 2 )/dQ 2 ) at exactly Q 2 = 
0. Since in elastic electron scattering the form fac- 
tor can be measured down to small but finite mo- 
mentum transfers Q 2 only one might think about a 
change of the form factor at very small Q 2 . Carlson 
[2Tj refitted a subset of the data of Bernauer [55] 
at very low Q 2 but finally got results in agreement 
with the radius from the complete fits if higher or- 
der terms in the form factor Q 2 expansion were 
included. Wu and Kao [23] tried a "thorn" at 
very small Q 2 equivalent to increasing again the 
3rd Zemach moment (?p)(2) so much that agree- 
ment from eq.Q results. However, here a danger- 
ous oversight occurs since eq.([T]) contains the first 
terms of the perturbative expansion series .24 only. 
Since the higher order moments diverge for such 
extreme modifications of the form factor at very 
small Q 2 cq.|l]) this expansion breaks down and is 
not applicable anymore. 

Recently the low Q 2 extrapolation has been formu- 
lated in the realm of chiral dynamics |25| putting 
it on firm theoretical grounds beyond the empirical 
Q 2 expansion. However, no change of the radius 
is indicated |26j . Similar ideas have been pursued 
also in ref. |27j . 

In a very recent paper Lorenz, Hammer, and 
Meifiner [28 try their own fit to the data of 
Bernauer et al. [5] [52] using the "continued frac- 
tion" model for the electric and magnetic form fac- 
tors Ge,m U- However, of the many models tried 
by Bernauer et al. this model was one of the worst 
delivering unstable fits due to poles outside the fit- 
ted Q 2 range. Only for Q 2 > 0.1 GeV 2 the fits of 
Lorenz, Hammer, and Meifiner are sufficiently sta- 



ble to allow an extrapolation to Q 2 = 0. However, 
the poles make such an extrapolation very ques- 
tionable and consequently also the radius derived 

m\- 

Also it has to be realized that the good agreement 
of the electronic Lamb shift and the electron scat- 
tering results make explanations with extrapolation 
modifications unlikely. 

• Dispersion relations 
In the already cited paper [28] Lorenz, Hammer, 
and Meifiner present also a refit of their older dis- 
persion relation fits complementing their old data 
base with the new data of Bernauer et al. [5] . As 
the old fits these new fits have a large X 2 e d = 2-2 
for 100 degrees of freedom (dof ) estimated from the 
figure; the dof are not given. Although the use of 
the theoretical x 2 distribution is not really justified 
for these fits since the theory function is unknown 
(Baysian situation) and the errors are not Gaussian 
distributed such a X 2 ed is equivalent to a probabil- 
ity P(Xr eti > 2.2; dof w 100) « 10~ 10 . The sensi- 
tivity of such fits to the radii is not discussed. 

• Fancy particles beyond the standard model 
In view of the convincing experimental evidence 
and the unsuccessful attempts with the mentioned 
conventional explanations some authors have spec- 
ulated about new physics. 

Batell, McKeen and Pospelov [3D] consider new vec- 
tor and scalar forces with carriers of less than 100 
MeV mass. Tucker-Smith and Yavin [3T] propose 
a new scalar or vector boson of about 1 MeV mass 
coupling more strongly to the muon than to the 
electron. Barger et al. [52] investigate the possibil- 
ity of new particles with a special coupling to the 
muon more generally, but conclude that new spin-0, 
spin-1 and spin-2 particles are disfavored by other 
experimental constraints. 

In a very recent paper Carlson and Rislow [33] dis- 
cuss two models, one involving new particles with 
scalar and pseudoscalar couplings, and a second in- 
volving new particles with vector and axial cou- 
plings. Though it appears not impossible to ac- 
commodate the new particles, some fine tuning of 
masses and couplings is needed to adjust them to 
the used constraints of the Lamb shift, muon mag- 
netic moments and kaon decay rate data. 

AN OVERLOOKED ISSUE 

The muonic Lamb shift is in leading order caused by 
the vacuum polarization in QED described by the ex- 
change of an electron-positron loop. The loop correc- 
tion can be calculated perturbatively with Feynman di- 
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agrams, see e.g. ref. [3U[3S], or non perturbatively from 
the exact solution of the Dirac equation using Green's 
function [36] [37] . Not surprisingly both approximations 
give the same result in leading order of the expansion 
in Za 2 . We stick to the first approach and note that 
it is equivalent to time dependent perturbation theory 
(tDEPpth). The energy shift is then given by the follow- 
ing expression: 

<Un< = -W;°I^M;o>' (3) 

|n/;0) is the unperturbed solution of the wave equation 
of the muon in the external Coulomb field. com- 
prises in principle the self energy, the vertex correction 
and the vacuum polarization propagator, but for muonic 
hydrogen it suffices to restrict oneself to to the exchange 
of one electron-positron loop, i.e. the vacuum polariza- 
tion. The formula is derived using Feynman rules and a 
first order expansion in energy with respect to the un- 
perturbed energy eigenvalue E^) ([31], chapter 14.2). 
Since the matrix element is of the form j T,* A p(r)r 2 dr 
with p(r) — \nl;0) 2 the charge distribution of the muon 
orbit, it is suggestive to define an external potential 
Vyp = — the "Uehling potential" . However, one has 
to be aware that this identification is not a trivial step 
but rather an approximation which has to be justified 
in the realm of QED. Weinberg shows how the external 
Coulomb potential emerges from the summation of multi 
photon exchange diagrams for a non relativistic two par- 
ticle system ([31], chapter 13.6). A similar derivation had 
to be performed here in order to show the validity of the 
external potential approximation. In this approximation 
the vacuum polarization propagator is used as the ex- 
ternal "Uehling potential" [38] . This point is essential 
here because it means one assumes that the muon can 
be regarded as a test charge in the "Uehling potential" 
[55] , The external "Uehling potential" is of order Za 2 
and does not include higher order contributions as two 
loop exchanges, i.e. "double vacuum polarization" [35] 
or "vacuum polarization iteration" [9 or "polarization 
insertion in two Coulomb lines" Carroll et al. [5J. 

Since tDEPpth breaks down for systems with station- 
ary bound states we have to solve a time independent 
wave equation ( 34 , chapter 14). This means one calcu- 
lates the effect of the external "Uehling potential" Vyp 
on the bound state eigenvalues and wave functions by 
solving the wave equation: 

H\nl) = E nl \nl) (4) 

with 

H = H + VyP = T + Vcoulomb + Vyp (5) 

where \nl) is now the new wave function and E n ; the new 
energy eigenvalue. If one solves eq.Q numerically this 
eigenvalue is exact or "non-perturbative" [BJ in the realm 



of time independent perturbation theory (tINDEPpth), 
but it stays first order in tDEPpth. Consequently, we 
shall call the energies and wave functions calculated with 
Vyp "perturbed" and without "unperturbed". We can 
easily bring eq.Q in a form corresponding to eq.([3]). Wc 
set \nl) = \nl;0) + S\nl) with (nl:0\nl;0) = 1 and 
5\nl) _L \nl;Q). Multiplying eq.Q with (nl;0\ fr om the 
left we get: 

E n i-E<$ ={nl;0\V V p\nl) (6) 

The normalization choice for \nl;0) and \nl) is the same 
as in conventional tINDEPpth yielding this simple form 
and making the situation lucid. However, one has to re- 
alize that the perturbed state \nl) is not normalized to 
1. As long as one does not make the relation to the un- 
perturbed state explicit this means just an overall factor 
leaving the energy eigenvalue unchanged. But here we 
do relate to E^ and one has to consider this carefully. 
(For a discussion of this point in the frame work of tIN- 
DEPpth (see [40], chapter 5.1.) 

We emphasize that the perturbed, i.e. exact, solu- 
tion contains only contributions of order Za (Coulomb 
potential) and Za 2 ("Uehling potential"). This is evi- 
dent if one looks at the expansion of these states in tIN- 
DEPpth [40]. The expansion coefficients are the matrix 
elements of the perturbation Vyp and mix the discrete 
unperturbed eigenstates. The idea that the electron- 
positron loop is reiterated in these solutions as in the 
Lippmann-Schwinger equation, i.e. tDEPpth, is not cor- 
rect for bound states [H], chapter 6.3. Stationary bound 
states cannot be reached in such an expansion since they 
are eigenstates with discrete eigenvalues. The point is 
that vacuum polarization loops occur only in QED and 
can be approximated in a wave equation only by an exter- 
nal potential, i.e. the "Uehling potential" here. There- 
fore, it is not correct to assume that the "polarization 
insertion in two Coulomb lines", Carroll et al. [6], is 
already present in the solution of the time independent 
wave equation eq. (|6| . Rather this contribution has to be 
calculated separately applying Feynman rules and added 
to the propagator E^. 

From this discussion follows that eq. ^ does not rep- 
resent a "non-perturbative" Lamb shift. In order to cal- 
culate the corrected Lamb shift with the perturbed wave 
function one has to propagate these wave functions with 
the propagator E^. In other words, one has to insert 
the approximate wave function gained form tINDEPpth 
with the external "Uehling potential" in the tDEPpth. 
One could think that this reinsertion is a double appli- 
cation of the Uehling effect, however, it rather realizes 
the dynamic effect on the energy shift due to the vac- 
uum polarization through the perturbed wave function. 
The forward-backward change between the "Uehling po- 
tential" as a propagator and an external potential is, of 
course, somewhat intuitive but we have no better ap- 
proximation. The procedure resembles in a way the self 
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consistent construction of the orbits of a multi electron 
atom. 

The derivation of this procedure is implied in the 
derivation of the Lamb shift in Weinberg's book ( |34j . 
chapter 14.2). Instead of expanding the perturbed solu- 
tions of the wave equation Un = un + Sun, or \nl) = 
\n I; 0) + S\n I) in the notation used here, with un the un- 
perturbed wave function and neglecting the terms with 
Sun one keeps the Un- Since the Uehling effect relates to 
the unperturbed eigenvalue E^j one has again to expand 
the energy denominator with respect to this energy. The 
result is a corrected formula for the Lamb shift -©Lamb,™ ; 
were one replaces the normalized unperturbed wave func- 
tion \nl;0) by the normalized perturbed \tiI)n in eq.(|3]): 

4amb,n! = N {n l\V V P \n I) N (7) 



provided by Pachucki 35J: 

V VP (r) = - — - ^ ^expi-m^utf) (8) 
where 




and q 2 is the internal momentum squared normalized to 
m 2 . For the radial unperturbed wave functions we use: 



This expression is intuitively plausible and has been al- 
ready mentioned by Carroll et al. [6], but, was not used 
and discussed by them. 

We can interpret eq. ([7| as the propagation of the per- 
turbed wave function \nl) by the vacuum polarization 
propagator, where, of course, the initial and final states 
have to be identical. If we compare eq.([6]) to eq.Q we 
realize that the perturbed state (n l\ appears only once in 
eq. ^. This means that one neglects the dynamical feed 
back effect of the perturbed wave function on the vacuum 
polarization, i.e. the additional energy shift caused by it. 

Carroll et al. [Bj take the difference E2 p — E2 s as 
the corrected Lamb shift with E2P and E2S the "non- 
perturbative" energy eigenvalues calculated by integrat- 
ing the Dirac equation numerically. At first glance one 
could believe that this is correct since the degenerate 



eigenvalues E^ — E 2 J = E 2 p cancel, but by comparing 
eq.Q and eq.([7]) we see that the dynamic feed back ef- 
fect is missing. We shall see in the next section that this 
explains the difference between the calculation of Carroll 
et al. [BJ and the evaluation of the Lamb shift with the 
perturbed wave functions. 



NUMERICAL REALIZATION 

In distinction to Carroll et al. in ref. [BJ we restrict our- 
selves to the non relativistic Schrodinger equation with 
the Coulomb potential of a point charge. This suffices for 
demonstrating the correction of the Lamb shift without 
considering fine structure splitting, finite size effect, etc.. 
Relativistic effects on the Lamb shift are small [9] and 
do not change anything essential for this discussion. This 
means we calculate the exact solutions of the Schrodinger 
equation with the point Coulomb potential and the ex- 
ternal "Uehling potential" added according to eq.([5]). 

A compact representation of the "Uehling potential" 
Vyp, well suited for our calculation, is the representation 



R2S = 



R 2 p = 



1 1 



1-1 
1 



2a 



2a ' 



1 r r 

— /= — ;= exp( ) — 

2^Q^4 K 2a J a 



(10) 
(11) 



with ao = he/ {an) the Bohr radius and \x the reduced 
mass. Using these wave functions in eq.([3]) we get for 
the leading order Lamb shift Ai?L am b = 



(o) 

Lamb. 2 P 



E 



(0) 



= 205.005 meV in agreement with 



Lamb, 2 S 

For the numerical integration of the Schrodinger equa- 
tion we have used Mathematica. As Carroll et al. [6] 
we have made extensive tests to guaranty the quality of 
the solutions. All calculations have been made with an 
internal precision of 64 digits and an accuracy goal of 
20 digits. The optimal method is the "Explicit Runge- 
Kutta" integration for the S-State and the change be- 
tween various methods provided "automatically" for the 
P-state. The numerical eigenvalues of the unperturbed 
2S and 2P states, i.e. without the "Uehling potential", 
are compared to the non relativistic exact solution (Bohr 
energies) and found to be good to a few neVs for differ- 
ent boundary conditions at small (~ 0.1 fm) and large 
(« 10000 fm) radii to which the eigenvalues are sensitive. 
Of course, since we do not take the difference of the large 
energy eigenvalues as Carroll et al. [B] but calculate the 
small Lamb shift only, we do not really need this extreme 
accuracy. However, when using our unperturbed numer- 
ical eigenvalues we get AELamb = 205.005 meV in com- 
plete agreement with the calculation using exact wave 
functions of eq.(10) and eq.(ll). The eigenvalues with 



the "Uehling potential" have been determined using the 
virial theorem. 

Figure [T] and Fig. [2] show the difference of the density 
of the unperturbed state \nl;0} 2 r 2 minus the density of 
the perturbed state \nl) 2 N r 2 . 



5 




1000 2000 3000 4000 5000 6000 
r[fm] 



FIG. 1. The difference (\2S;0) 2 - \2S) 2 N )r 2 showing the po- 
larization charge density divided by the negative elementary 
charge of the muon due to the "Uehling potential" for the 2 S 
state 
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FIG. 2. The difference (|2P;0) 2 - |2 P)%)r 2 showing the po- 
larization charge density divided by the negative elementary 
charge of the muon due to the " Uehling potential" for the 2 P 
state 



One nicely sees the polarization charge induced in the 
vacuum by the muon. The wiggles at small radii in Fig. 
[I] are no numerical artifacts but due to the double bump 
structure of the 2S state. In order to get an idea of 
the scales we note that the Bohr radius for the muon is 
ao = 285 fm, the scale of the Compton wave length of the 
electron or positron in the electron-positron pair of the 
vacuum polarization A = 2426 fm, the rms radius of the 
2S state (2 S\r 2 \2 S) 1 / 2 = 1854 fm, and for the 2P state 
(2P|r 2 |2P) 1 / 2 = 1560 fm. As expected some positive 
charge is pushed to larger radii compensated by negative 
charge at small radii indicating the induction of the po- 
larization cloud in the vacuum. An analytical derivation 
of this polarization is derived in ref. [51], chapter 11.2. 

Calculating the difference of the Lamb shifts with the 
normalized perturbed 2S and 2P states according to 
eq.([7]) one gets the salient result of this paper: 



A ££amb CharS ° = 205.307(l)meV 



(12) 



the value with different integration boundaries. Compar- 
ing this result to the canonical value for the unperturbed 
wave functions AE^^ lnt chalgc = 205.005(1) one arrives 
at 



5(AE hamh ) =0.302(1) meV 



(13) 



in very good agreement with the searched for difference 
of eq. (U). 

For the difference of the eigenvalues from the numeri- 
cal solution of the wave equation according to eq.([6]) we 
get E 2P - E 2S = 205.156(1) meV. If we take directly 
the eigenvalues E 2 p — E 2 s from the numerical solution 
of the wave equation eq.Q as Carroll et al. [5] we get 
205.159(3) meV. The deviation from the difference calcu- 
lated with eq. ^ indicates a limit of the numerical accu- 
racy. Since we do not use the difference of the eigenvalues 
we have not insisted to improve this limit. Considering 
the relativistic correction due to the Dirac wave functions 
of 0.021 meV |35j . missing in our non relativistic calcula- 
tion, this is in good agreement with the relativistic result 
of Carroll et al. |6j of 205.1706(5) meV. 

As already discussed this result does not include the 
higher order Za A contribution as the two loop vacuum 
polarization exchange. Therefore, the "double vacuum 
polarization" contribution of 0.151 meV [35] has to be 



added to our result (5(A£ , Lamb) in eq. (13) as well as to 
that of Carroll et al. [5] of 205.1706(5) meV, who, how- 
ever, neglected it. 



CONCLUSIONS 

If the Lamb shift is taken as the dynamical QED ef- 
fect due to the interaction of the muon with the vacuum 
polarization and not as a shift caused by the "Uehling 
potential" as an static external potential, one gets agree- 
ment for the radius determined from the Lamb shift in 
muonic hydrogen with the combined electronic experi- 
ments. Since the relativistic calculations including the 
finite size effects of Carroll et al. [B] have to be redone 
realizing the new considerations of this paper we stick to 
the formula eq. ([T]) used by Pohl et al. [T]. This means 
all other corrections, in particular the "double vacuum 
polarization" and the finite size effect, are the same as 
used in that analysis. If we correct the Lamb shift of 
the point charge in eq.|l]) to the value calculated in this 
paper, we arrive at a new value for the rms radius of the 
proton derived from the muon experiment: 



(r 



2\l/2 



0.87650(71) fm 



(14) 



where we have taken the recent 3rd Zcmach moment from 
ref. [16] and included the estimated error of the point 
charge Lamb shift. This value is now in good agreement 
with the best electron scattering rms radius 0.879(8) fm 



where the error is a best estimate from the variation of 



Again QED wins and there is no reason to fear a " chink 
in the armor" as in ref. [42]. However, it is also true 
that "it must be repeated that the theory of relativistic 
effects and radiative corrections in bound sates is not yet 
in entirely satisfactory shape" [33], p. 560. 

I am indebted to my son Johannes Walcher for hints 
and remarks and to Marc Vanderhaeghen and Vladimir 
Pascalutsa for a critical discussion of the prefmal paper. 
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